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We show that the breaking time of quantum-classical correspondence depends on the type of 
kinetics and the dominant origin of stickiness. For sticky dynamics of quantum kicked rotor, when 
the hierarchical set of islands corresponds to the accelerator mode, we demonstrate by simulation 
that the breaking time scales as r% ~ with the transport exponent ji > 1 that corresponds 

to superdiffusive dynamics JlTj . We discuss also other possibilities for the breaking time scaling and 
transition to the logarithmic one Tn ~ ln(l/7i) with respect to Ti. 
PACS number: 05.45Mt 

1. Classical chaotic dynamics can be characterized by a Lyapunov exponent A and infinitely divisible filamentation 
of the phase flux. Quantized procedure stops the filamentation due to the uncertainty principle and, as a result, 
breaks the applicability of semiclassical approximaiton. The corresponding breaking time was found in § 

7fc = (l/A)In(J /n) ) (1) 

where Iq is a characteristic action, indicating a fast (exponential) growth of quantum corrections to the classical 
dynamics due to chaos. The origin of this time was explained in detail in Q and gained a wide discussion ||. The 
logarithmic scaling in h for corresponds to a fairly good and uniform chaotic mixing (see also a final comment in 

ft- ... . 

Typical Hamiltonian chaotic dynamics is not ergodic due to the presence of infinite islands set in phase space [g[ , 

and Lyapunov exponent is not uniform due to cantori and possible hierarchical structures of islands ||] and their 
stickiness [^J. This complicates a process of diffusion transforming the transport from the Gaussian type to the 
anomalous (fractional) one ]8|,^| . Particularly, sticky properties of the island boundaries should impose algebraic laws 
of the survival probability 

P(t)~lA 7 (2) 

of a particle to escape after time t from a domain near the islands HQ . The immediate consequence from the scaling 
property (|^) is its breaking at some critical time t* for the case of quantum chaos since a hierarchical dynamical 
chain has no limit t — > oo and should be abrupted when an island in the chain reaches area S* = fi and the quantum 
effects become important. This comment, starting from JTo[ , was discussed in detail in p|, where a power law for the 
Planck constant scaling was suggested 

t>: ~ i/n s (3) 

for the breaking time of classical considertion applicability with a value of 5 depending on the type of the classical 
algebraic escape time distribution. There were no quantum simulation in these works and no definite values of i5. 
Discussion of the algebraic law (^) started in [§,0, was continued in recent publications pd|-Jl6[| on the basis of 
simulations of quantum maps with an explicit evaluation of (^J). 

The questions arise: What is the actual scaling of breaking time of quantum-classical correspondence with respect 
to the Planck constant, logarithmic or algebraic? How universal are values of 7 or 81 Using some results of |l7],[l8| on 
the strong derealization effects, and simulation for the quantum kicked rotor (QKR) we will show here that the result 
(||), being valid for the cases of normal (Gaussian type) diffusion, appears to be an algebraic of the type @ when 
diffusion becomes anomalous, i.e. superdiffusion, with the second moment of the truncated distribution function as 

<P 2 ) ~ f (4) 

and transport exponent fj, > 1. We were able to obtain t£ from the simulation as a crossover time of the survival 
probability that changes the exponent of its algebraic behavior, and compare values of 5, predicted by the theory 
with the one obtained from the simulation. 

2. We consider QKR that corresponds to the standard map in the classical limit 

p' = p + K sing , q' = q + p' (5) 
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defined on the cylinder p S (—00, 00), q S {—if, 7r) with a control parameter if and the Lyapunov exponent A ~ 1/ h\K 
for K ^> 1 and for almost all domain excluding areas where 2f| cosq| < 1. The marginally stable points are defined 
by the conditions K m = 2irm, p = 2irn, q — ±ir/2 with integers (m, n). For 

< K - K m < AK m (6) 

a new set of islands appears f6] Jl9| , pof , called tangle islands in |2(|, as a result of bifurcation. Dynamics inside the 
islands is known as the accelerator mode |2l],|2^] and we will call them accelerator mode islands (AMI). Changing 
of K within the interval (@) influences strongly the topological structure of AMI and consequently the values of the 
transport exponent fi — (jl[K) [^3|Js[|, since the stickiness of trajectories to the islands boundaries can be different for 
different island topologies. 

It was established in || that for a special ("magic") value of K = K* = 6.908745. . . the stickiness can be especially 
what makes this case to be convenient to study properties of the anomalous transport. For the K* it appears a 
hierarchical set of islands- around-islands with the islands sequence 3 — 8 — 8 — 8 — .... The islands chain satisfies the 
renormalization conditions 

S^+V = A s 5 (n) , T(™ +1 > = A T T (n) , 

jV^ 1 ) = \ N N^ , (7) 

where n is a number in the hierarchy sequence, S^ n ' is an island area, 

T (n) 

is a period of the last invariant curve of the 
corresponding island, N^ n ' is a number of islands in the chain of the n-th hierarchy level, and As < 1, Ax > 1, Ajv > 1 
are some scaling parameters. The renormalization transform (0) can be also extended for Lyapunov exponents A in 
a sticky area of the island's boundary of the n-th level of the island's hiearchy p3j : 

A(" +1) = AzAW = a^AW , (\ L < 1) • (8) 

In the absence of the island's hierarchy, we get just the result ([j]) with A = A^ ^ . In the presence of the island's 
hierarchy we can introduce particle flux in phase space through the island's chain of the n-th hierarchical level. 
It reads 

$(n) _ s^N^/T^ = $(°) (As Aw/At)" , 

$(°) = 5(o)jv(°)/tW (9) 

in correspondence to (Q) and (|J). For K = K* and the corresponding island's hierarchy A^r = At [[§ and thus 

*W = Ag$W . (10) 

The quantum mechanical consideration of the proliferation of islands is meaningful until the smallest island size is 
of the order of fi. Therefore we get from (0) , @ , ( 10|) : 

q . — f, — q( n «) — q( a )\ n o 
$ min = $(°)AS° = A" S* (0) Ar(°)/T( ) = hNW/TW , (11) 
and the quantum "cut-off" of the hierarchy appears at 

n = |ln(?i/5( ))|/|lnAs| = |ln/ l |/|lnAs| , (12) 

where we introduce a dimensionless semi-classical parameter h = Ti/Iq and specify S^ 1 = la. Particularly, for the 
hierarchy at K — K* we have — 3 and A at = 8 but it could be many other hierarchies (see more in |2J]). After 
the substitution A = A( n °) we get with (§) and (0): 



r n = (1/A("))) \n(l/h) = (1/h) 1 ^ \n(l/h) (13) 

with 

H= |lnAs|/lnA T ■ (14) 

The expression of \i through As, At was found in for the considered islands hierarchy and the expression for r% 
coincides with the obtained in |L7| in a different way. The expression (|l3|) is close to (|^) up to a logarithmic term and 
defines 
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5 = 1/u, (15) 

which for K = K* provides \x = 1.25 (see Q) and d = 0.8. When there is no islands hierarchy, we may put As — > or 
uq — > 0. This yields transition from the algebraic law ( |l3| ) to the logarithmic one for the breaking time and resolves 
the paradox discussed in Q. In addition to this, it was shown in Q that for the considered islands hierarchy 

7 = l + /i = l + |lnAs|/lnA T , (16) 

what gives 7 = 2.25 for K = K*. Just these values of 7 and 8 we have checked by a simulation. 

3. Numerical study of the problem is based on investigation of the quantum survival probability in some domain 
Ap G (— 7r,7r) that includes the islands hierarchy f25|| . The main difficulty appears due to a part of the wave function 
that belongs to an island interior and "flies" fast along p. A simple way to avoid this type of the "escape" from Ap 
is to apply a shift operator 

J = cxp{-(27r /h)d/dh} = exp{-2mq/h) , (17) 

where q is a dimensionless coordinate, p — hn — —ihd/dq is a dimensionless momentum operator and the wave 
function ty t at a discrete time t is considered in the coordinate space. In analogy to jll];^6) we also introduce the 
absorbing boundary conditions at the edges of the interval An € [—(N+l)/2, (N—l)/2] for the momentum eigenvalues 
hn. Then the quantum map that keeps the information of the trapping into An part of the wave function is 

9 t+ i=-pjU*t (18) 

with the evolution operator 

U = exp ( - ihh 2 /2) exp ( - i(K/h) cos q) . (19) 

As usual, for the numerical convenience the dimensionless Planck constant h is taken in a form h = 2it/(N + g), where 
g= (V5- l)/2 is the inverse golden mean. 
The survival probability is defined as 

P(t) = |* t | 2 (20) 

together with definition (|T8|). Simulation was performed for 12 values of N from the interval (5 • 10 3 -j- 7.5 • 10 4 ) which 
corresponds to a good semi-classical approximation for a fairly long time until it fails. 

A typical behavior of P(t), obtained for K — K*, is shown in Fig. 1. It consists of the crossover from the classical 
behavior @ with 7 w 2.25 (the same as in ||) to some very different dependence. The crossover points Tn can be 
identified in a way as it is shown in Fig. 1 for different values of N, i.e. h. The corresponding result is presented in 
Fig. 2. It provides the value of 1/5 w 1.33 ± 0.36 w \x in a good agreement to the presented theoretical estimation 
(see©). 

4. In conclusion, we need to underline that there is no unique scenario of chaotic diffusion in classical limit and 
therefore one may expect no unique breaking mechanism of quantum classical correspondence. Stickiness and algebraic 
kinetics through cantori was considered in [pi with a specific choice of the Markov tree that defined a type of kinetics 
and corresponding scales. In this consideration a special hierarchical set of resonance islands was selected, while in 
P,p3| the Markov tree was constructed for the tangle islands (see more about the classification in p(J ) . The difference 
in the choice of the islands set selection is imposed by the value of K. In this article we choose K — K* which leads 
to dT|) while in |l|,||,|| the value K < 2n was selected which probably leads to stickiness phenomenon described 
in (§] and to the value 

5=1/7. (21) 

This value was not linked to the transport exponents u. Let us mention also the value 5 <~ 0.5 proposed in p6[ for 
K = 2.5 when sticky islands set appears without accelerator mode and with no superdiffusion. 
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FIGURE CAPTIONS 

Fig. 1 Typical evolution of the quantum survival probability for N = 25557. The dashed lines, which correspond to 
the asymptotics slopes, determine the breaking point r^. 

Fig. 2 The quantum breaking points t% vs dimensionless semiclassical parameter h. The solid line with the slop 
1/8 « 1.33 ± 0.36 corresponds to least square calculations, and the dashed line is the analytical prediction with 

6 = 
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